Centrality dependence of ipt over J/ip : a signal of QGP 
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If nuclear collisions lead to QGP formation then the ratio 
of ipt over J/ip will remain constant with Et as both J/ip 
and t/i/ are melted in the QGP. On the other hand, if hot 
hadronic matter is produced, the ratio will continually fall 
with Et, as ipt's are more suppressed in hadronic matter than 
the J/ip's. We have constructed the ratio for Pb+Pb collisions 
at SPS from the existing NA50 data. The ratio gives the 
indication of a possible QGP formation at SPS energy, but 
definite conclusion can not be reached. We have also given 
the prediction for the ratio at RHIC energy. 



Lattice QCD predicts that under certain conditions 
(sufficiently high energy density and temperature) ordi- 
nary hadronic matter (where quarks and gluons are con- 
fined) can undergo a phase transition to deconfined mat- 
ter, commonly known as quark gluon plasma (QGP). J/ip 
suppression is recognized as one of the promising signal of 
the deconfinement transition. Due to screening of color 
force, binding of cc pairs into a J/ip meson will be hin- 
dered, leading to the so called J/ip suppression in heavy 
ion collisions [1]. Experimental data indeed show sup- 
pression. However, all the data prior to NA50 Pb+Pb 
are well explained in terms of nuclear absorption, also 
present in pA collisions [2,3]. NA50 collaboration [4] ob- 
served anomalous suppression (i.e. suppression beyond 
the normal nuclear absorption) in 158 GeV/c Pb+Pb 
collisions [4]. The ratio of J/ip yield to that of Drell- 
Yan pairs decreases faster with Et in the most central 
collisions than in the less central ones. It has been sug- 
gested that the resulting pattern can be understood in a 
deconfinement scenario in terms of successive melting of 
charmonium bound states [4]. However, it was realized 
later that the data could be well explained in a variety of 
models [5-9], with or without QGP formation. In ref [8] 
it was also shown that the predicted J/ip suppression at 
RHIC in a QCD based nuclear absorption model, agree 
well with predictions obtained in a QGP based model. It 
seems that, even at RHIC, deconfining phase transition 
could not be detected from the J/ip suppression. 

Not only J/tp, but other states of charmoniums (e.g. 
X and ipt) are also suppressed in a QGP or in a nu- 
clear matter. In pA collisions, ipt suppression is simi- 
lar to J/ip suppression [10]. Recently in Quark matter 
2002, NA50 collaboration confirmed that in pA collisions, 

a abs' N ~ a abs 1^1 ■ They measured J/ip as well ipt from 
Be, Al, Cu and W targets. Parameterizing the produc- 



tion cross section as a pA cx A a , NA50 collaboration ob- 
tained, a J/v , = 0.933 ± 0.105 and ay,, = 0.906 ± 0.022. 
Nearly identical values of a for both the J/ip and ipt, 
in a Glauber type of model of nuclear absorption trans- 
late into similar absorption cross section for them, i.e. 
a abs' N ~ °abs ' ' contrary to the popular expectation that 
ipt being twice as large in size than the J/ip, cC 6 ' s will be 

much larger than cr Q Y<f W - The apparent contradiction is 
resolved in the color octet model [12]. In the color octet 
model, perturbatively produced cc pairs first neutralizes 
its color by combining with a soft collinear gluon. The 
pre-resonance ccg state then transforms in to a proper 
charmonium state, J/ip or ipt. In pA collisions, the nu- 
clear medium sees only the prc-rcsonancc state. Equality 
of &alf N and cr^/// is then explained. 
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Unlike in pA collisions, in AA collisions, J/ip and 
ipt suppression differs. NA38/NA50 collaboration mea- 
sured centrality dependence of J/ip as well as of ipf's in 
S+U/Pb+Pb collisions [13,4]. Data indicate that com- 
pared to J/ip, ipt's are more suppressed. For example, 
in S+U collisions, from peripheral to central collisions 
J/ip's are suppressed by a factor of ~ 1.3, while the ipt's 
are suppressed by a factor of ~ 4. Similarly in Pb+Pb 
collisions, while ipt's are suppressed by a factor of 8, J/tp's 
are suppressed by a factor of three only. Thus in AA col- 
lisions, additional suppression mechanism is operative for 
ipt's, which is absent for the J/ip's in AA collisions. 

One of the source for additional suppression could be 
the QGP formation. If QGP formation is the source of 
the additional suppression of ipt's, why the effect is not 
seen in J/ip? Color screening studies shows that in a 
QGP, melting of J/ip require a temperature of 1.2T C , 
while the ipt's are melted at T c only. Thus if QGP is 
produced around T c , its effect may be felt only on ipt, 
not on J/ip's. Also as the time scale of production of 
ipt is less than that of J/ip, ipt's can better probe the 
initial condition of the produced matter. Thus it is pos- 
sible that effect of QGP formation will be seen only in 
ipt rather than in J/ip. Hadronic comover's could also 
be the source of additional suppression. In AA collisions 
a large number of secondaries are produced. Absorption 
cross section of ipt in comovers could be larger (due to 
larger radius) than that of J/ip's, leading to increased 
suppression of ipt. 

In the present letter we have analyzed the NA38/NA50 
data on the centrality dependence of ipt over Drell-Yan 
ratio, in S+U and in Pb+Pb collisions. Analysis shows 
that absorption in comovers or in QGP, both the sce- 
nario could explain the data. Even at RHIC energy, the 
ambiguity is not removed. It may not be possible to 
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detect deconfinement phase transition from J/ip or ip/ 
suppression. Next we consider centrality dependence of 
the ratio of ip/ over J/ip. Gupta and Satz [14] proposed it 
as a signal of QGP. The idea is simple. If in the collision, 
QGP is formed above a temperature 1.2T C , then both 
the J/ip and ip/ will be equally suppressed and the ratio 
will remain constant with Et- Otherwise, the ratio will 
continually decrease with Et, as ip/ are more suppressed 
than J/ip. The ratio has been considered as a thermome- 
ter for the deconfinement temperature also [15]. Though 
simple and quite old idea, unfortunately, NA50 collab- 
oration did not present their results for the said ratio 
for Pb+Pb collisions, which generated so much interest 
about possible deconfinement phase transition. From the 
existing data, we have constructed the ratio for Pb+Pb 
collisions. The centrality dependence of the ratio, though 
show a tendency towards saturation, it is not possible to 
conclude decisively about phase transition. Conclusive 
signal could be obtained at RHIC energy. 

In the QCD based nuclear absorption model [8,16], 
J/ip production is assumed to be a two step process, 
(a) formation of a cc pair, which is accurately calcula- 
ble in QCD and (b) formation of a J/ip meson from the 
cc pair, which is conveniently parameterized. The J/ip 
cross section in AB collisions, at center of mass energy 
,/s is written as, 

a J ^{s)=KY J J dq 2 (%^) / dx F <p a/A (x a ,Q 2 ) (1) 

a,b 

/B^Q^^^xF^j^q 2 ), 



q 2 + 



where J2 a b runs over an P ar ton flavors, and Q 
Am 2 . The K factor takes into account the higher order 
corrections. The incoming parton mome ntum fractio ns 
are fixed by kinematics and are x a = {\J x 2 F + AQ 2 /s + 
x F )/2 and x b = (y/ 'x 2 F + AQ 2 /s - x F )/2. a ab ^ cS 
are the subprocess cross section and are given in [17]. 
F C c-*j/4)(q 2 ) is the transition probability that a cc pair 
with relative momentum square q 2 evolve into a physical 
J/ip meson. It is parameterized as, 



relative square momentum gain per unit length. Param- 
eters of the model (ap,KNj/^ and e 2 ) can be fixed from 
experimental data on total J/ip cross section in pA/AA 
collisions. In Fig.l, NA50 high statistics data [11] are 
shown. Both the data sets are well explained in the model 
with £ 2 =0.1875 GeV 2 / fm. Nuclear suppression of J/ip 
and ipi are due to same mechanism, i.e. gain in the rela- 
tive 4-square momentum of the cc pairs. Naturally, J/ip 
and of ip/ shows similar A-dependence. 
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FIG. 1. (a) The experimental ratio of total J/tp cross sec- 
tion and Drell-Yan cross sections in pp and pA collisions. The 
solid line is the fit obtained in the QCD based nuclear absorp- 
tion model, (b) same as (a) for ipt 

In our earlier work [18], we have shown that the model 
could reproduce centrality dependence of J / ip over Drell- 
Yan ratio in S+U and Pb+Pb collisions. The J/ip or ipi 
cross sections at an impact parameter b as a function of 
Et can be written as, 



dE T <Pb 



NN 



J d 2 sT A (s)T B (b-s)S(L(b,s))P(b,E T ) 

(3) 



where Ta.b is the nuclear thickness function, T A ,B{b) = 
J dzp(b,z). For the density, we have used the three pa- 
rameter Fermi distribution [19], 

po(l +^?4) 



F c ^ J/ 4q 2 )^N J/ ^q 2 )e(Am' 2 -Am 2 c -q 2 ) (2) P(0 = j + Z P {-{r^C)/a) ' j = A ( 4 ) 



(1 



4m' 
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In a nucleon- nucleus/nucleus- nucleus collision, the pro- 
duced cc pairs interact with nuclear medium before they 
exit. It is argued [16] that the interaction of a cc pair with 
nuclear environment increases the square of the relative 
momentum between the cc pair. As a result, some of 
the cc pairs can gain enough relative square momentum 
to cross the threshold to become an open charm meson. 
Consequently, the cross section for J / ip production is re- 
duced in comparison with nucleon-nucleon cross section. 
If the J/ip meson travel a distance L, q 2 in the transition 
probability is replaced to q 2 — > q 2 + e 2 L, e 2 being the 



The parameters of the density distribution (C,u) and 
a) are taken from [19]. 

P(b, Et) is the Et — b correlation function. We have 
used the Gaussian form for the Et — b correlation, 

P(b, E T ) cx exp(-{E T - qN p {b)f /2q 2 aN p (b)) (5) 

where N p (b) is the number of participant nuclcons at 
impact parameter b. a and q are parameters related 
to dispersion and average transverse energy. For S+U 
collisions at 200 GcV/c, the parameters are, o=3.2 and 
g=0.74 GeV [2], and for 158 GeV/c Pb+Pb collisions the 
parameters are, a=1.27 and g=0.274 GeV [5]. 
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In Eq.3 S(L) is the suppression factor due to passage 
through a length L in nuclear environment. At an impact 
parameter b and at point s, the transverse density is 
calculated as, 

n(b, s) = T A (s)[l - e -^NNT B (b- s )^ + [ A ^ B } ( 6 ) 

and the length L(b, s) that the J/ip or ipl meson will 
traverse is obtained as, 



L(b,s) = n(b, s)/2p 



(7) 



The Drell-Yan pairs do not suffer final state interac- 
tions and the cross section at an impact parameter b as 
a function of Et could be written as, 



2„DY 



So 



dE T d 2 b 



agl J d 2 sT A ( S )T B (h- S )P(b,E T ) 



(8) 



The additional suppression required for ipl may be due 
to QGP formation or due to comover interactions. To 
take into account the suppression due to QGP formation 
we assume that above a threshold density, n c , all the ipl 
are dissolved [5], and introduce the additional suppres- 
sion factor Sqcp in Eq.3, 



<SqgpO,s) = 6(n c 



< E T > (b) 



i(b,s)), (9) 



The additional suppression factor in the comover sce- 
nario can be written as [20] , 

S co (b,s) = exp(-<7 co v re in (h,s)T \xi(R T / v re iT )) (10) 

In the above equation, a co is the comover absorption 
cross section for the ip/'s, v rei =0.6 is the relative veloc- 
ity of ipl with respect to comovers and To =2 fm, is the 
time beyond which the comover interactions starts. Rt 
is the transverse radius of the system and n is the initial 
comover density. To account for the variation of density 
with E T , we take n =< n > E T / < E T > (b = 0), 
with < n >=0.8 fm~ 3 [20]. The only quantity to be 
fixed is the a COl which we obtained directly from fitting 
the S+U data. 

The other unknown quantity is the N J ^^' — 
B w a( i{DY)'ll N ■ Experimentally it is known for 450 GeV 
pp collisions [11]. Craigie parameterization [21] of DY 
cross sections could be used to obtain its value at other 
energies. For 200 GeV/c S+U collisions, the extrapolated 
values are, Y J /'''=32-42 and A^'=0.53-0.68, for the DY 
invariant mass in the ranges of 2.1-3.1 GeV. We obtain 
the values of N J ^ and N^' from a constraint fit to the 
NA38 data such that they are within the range of extrap- 
olated values. For J/ip's in Pb+Pb collisions, we rescale 
the value by the factor 1.051 [11]. For ipl's, we use the 
same value. 

In Fig. 2a, we have shown the NA38 data on the cen- 
trality dependence of J/ip over Drell-Yan ratio, for 200 
GeV/c S+U collisions. The solid line is the fit to the 



data obtained with B lx ^a J J^ [ /<rj^=39. 02. Data are well 
explained. In Fig. 2b, the latest NA50 data on the cen- 
trality dependence of the ratio of J/ip over Drell-Yan are 
shown. The solid line is the prediction in the QCD based 
nuclear absorption model, with the normalising factor, 
B^a J J^/a^ Y N =A1.0l. The latest NA50 data are also 
well explained in the model. We note that there is no 
scope for additional suppression due to comover interac- 
tion or due to QGP formation. Thus centrality depen- 
dence of J/ip suppression in S+U or in Pb+Pb collisions 
do not require additional suppression due to QGP for- 
mation or due to comover interactions. 
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FIG. 2. (a) Experimental data on the centrality depen- 
dence of J/ip over Drell-Yan ratio, in 200 GeV/c S+U col- 
lisions. The solid line is the fit obtained to the data in the 
QCD based nuclear absorption model, (b) same as (a) for 
Pb+Pb collisions. (c)Experimental data on the centrality de- 
pendence of ipl over Drell-Yan ratio in S+U collisions. The 
solid, dashed and dash-dotted lines are the fit to the data with 
nuclear, nuclear+comover and nuclear+QGP suppression re- 
spectively, (d) same as (c) for Pb+Pb. 

The centrality dependence of ipl over Drell-Yan ra- 
tio on the other hand require additional suppression. 
In Fig. 2c and 2d the NA38/NA50 data on the ratio 
B^a(ipl)/a(DY) for 200 GeV/c S-U and for 158 GeV/c 
Pb+Pb collisions are shown. The solid line is the ratio 
obtained in the QCD based nuclear absorption model. 

Irrespective of £?uu— frr value, the model clearly fails to 

explain both the data sets, ipl's are not sufficiently sup- 
pressed to agree with experiment. As discussed in the 
beginning, additional suppression could be either due to 
comovers or due to QGP formation. In Fig. 2c and d, the 
dashed line is the ratio obtained with nuclear+comover 
suppression, with cr co =8 mb. For both the data sets, we 

have used S„„^f =0.59, obtained from fitting the NA38 

S+U data. The comover scenario fits the Et dependence 
of ipi in S+U and in Pb+Pb collisions. For Pb+Pb colli- 
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sions, for the very peripheral collisions, model produces 
more suppression than in data. Since we have fixed the 
comover density in central collisions, the simple ansatz 
may be inaccurate for peripheral collisions. 

In S+U collisions, centrality dependence of ipl over 
Drell-Yan ratio is not well explained if the nuclear sup- 
pression is augmented with suppression due to QGP for- 
mation (the dash-dotted line). At low Et data are not 
explained. Also we obtain a threshold density, n c =1.8 
fm~ 2 , which is too low for QGP formation. For Pb+Pb 
collisions on the other hand (Fig. 2d), centrality depen- 
dence of ipl over Drell-Yan ratio are rather well explained 
with nuclear plus QGP suppression. In Fig. 2d, the dash 
dotted line is the ratio obtained with threshold density, 
n c =2.8 fm~ 2 . Data are well explained throughout the 
Et range. However, as nuclear plus comover suppres- 
sion also explain the data, it is not possible to conclude 
positively about the formation of QGP from the Et de- 
pendence of ipl suppression. 

With RHIC being operational, it is interesting to pre- 
dict suppression at RHIC energy. At RHIC energy, the so 
called hard component, which is proportional to the num- 
ber of binary collisions, appear. Model dependent calcu- 
lations indicate that the hard component grows from 22% 
to 37% as the energy changes from 56 GeV to 130 GeV 
[22]. In our calculation, we have used 37% hard scatter- 
ing component. 

In Fig. 3, we have shown the predicted centrality de- 
pendence of the ipl over Drell-Yan ratio at RHIC energy 
for Au+Au collisions. The solid and dashed lines corre- 
sponds to nuclear+comover and nuclear+QGP absorp- 
tion respectively. They agree closely with each other. 
Et dependence of ipl over Drell-Yan ratio at RHIC also 
could not distinguish between the two scenarios. 
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FIG. 3. Centrality dependence of ipl over Drell-Yan ratio 
at RHIC. The solid and dashed lines are obtained with nu- 
clear+comover and nuclear+QGP suppression respectively. 

Next we consider the centrality dependence of ipl over 
J/ip. As told earlier, it has been proposed as a signal 
of the QGP formation. The proposal follows from the 
simple observation that in a QGP both the J/ip and ipt 
will be melted. Consequently, the ratio a (ipl) /a( J/tp) 
will remain constant with Et- Otherwise, the ratio will 
continually fall with Et, as ipl are more suppressed than 



J/ip. In Fig. 4, we have tested the proposition. NA50 
collaboration did not present the data. From the exist- 
ing J/ip and ipl data we have constructed the ratio. It 
is shown in Fig. 4. The ratio, for S+U collisions is also 
shown in Fig. 4. For S+U collisions, the ratio fall contin- 
uously with Et- QGP is not formed in the collisions. For 
the Pb+Pb collisions, the ratio falls with Et till around 
70 GeV and thereafter shows a tendency of saturation. 
Data do not cover enough Et range for a definite con- 
clusion. In Fig. 4, the solid and the dashed lines are the 
ratio for Pb+Pb collisions in the nuclear+comover and 
nuclear+QGP suppression. As expected both of them 
fits the data. We note that even at large Et, difference 
between the two model calculations is small (1-2%). Even 
if there is a phase transition, it will be difficult to reach 
a definite conclusion. 

The situation is much better at RHIC energy. Our pre- 
diction for the ratio at RHIC is shown in Fig. 4. The dash- 
dot line is the prediction obtained with nuclear+comover 
absorption. It shows continual fall of the ratio. In con- 
trast, with nuclear+QGP suppression (the dash-dot-dot 
line), the ratio remain constant for Et > 70 GeV. The 
difference between the two predictions is also large and 
easily detectable. 
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FIG. 4. Et dependence of the ratio, ipl over J/tp in 200 
GeV/c S+U and in 158 GeV/c Pb+Pb collisions. The solid 
and dashed lines are obtained with nuclear+comover and nu- 
clear+QGP suppression. The predicted ratio at RHIC energy 
with nuclear+comover suppression is shown as the dash-dot 
line. The dash-dot-dot line is the ratio with nuclear+QGP 
suppression. 

To conclude, we have analyzed the centrality depen- 
dence of J/ip and ipl suppression in S+U and in Pb+Pb 
collisions. It was shown that while the J/tp suppression 
is well explained in the QCD based nuclear absorption 
model, the model could not explain the centrality depen- 
dence of ipl suppression, ip/'s require additional suppres- 
sion, either due to QGP formation or due to comovers, 
two scenarios could not be distinguished, even at RHIC. 
We then considered the Et dependence of the ratio of ipl 
over J/ip as a signal for the deconfining phase transition. 
If QGP is formed following a deconfinement phase tran- 
sition, the ratio will remain constant with Et in contrast 
to the continuous fall of the ratio in case of no such for- 
mation. The experimental a(ipr)/a(J/ip) in Pb+Pb col- 
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lisions is not conclusive. However, at RHIC energy, the 
ratio could distinguish between the comover and QGP 
suppression. 
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